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been  proposed  as  discriainants  between  earthquakes  and  underground  nuclear  ex- 
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ABSTRACT 

Short-period  P-waves  first  arrivals  recorded  at  the  Large  Aperture  Seis- 
mic Array  in  Montana  from  a large  number  of  teleseismic  events  were  filtered 
in  the  bands  of  O.A-3.0,  0.4-0. 8,  1.4-1. 8,  1. 8-4.0  Hz.  Estimates  of  signal 
amplitude  in  each  band  were  used  to  form  spectral  ratios  similar  to  those 
which  have  been  proposed  as  discriminants  between  earthquakes  and  underground 
nuclear  explosions.  Values  for  the  ratios  scattered  over  nearly  two  orders 
of  magnitude  and  apparently  did  not  depend  on  earthquake  size  or  epicentral 
distance.  Average  logarithms  of  spectral  ratios  for  166  separate  geographic 
regions  ranged  over  one  magnitude  unit  and  showed  some  correlation  with  known 
patterns  of  lateral  Q inhomogeneity  in  the  earth,  suggesting  that  attenuation 
causes  much  of  the  spectral  ratio  scatter.  However,  attempts  to  correct  LASA 
m^  on  the  basis  of  the  observed  spectral  ratios  failed  to  reduce  the  scatter  in 
an  plot.  Variation  in  stress  drops  and  other  physical  processes  connected 

with  the  mechanism  of  earthquakes  must  cause  as  much  variability  in  observed 
spectral  discriminants  as  variations  in  attenuation.  Spectral  ratios  for  sev- 
eral explosions  in  the  USSR  did  not  separate  well  from  those  of  the  global  pop- 
ulation of  earthquakes  and,  more  Importantly,  they  did  not  separate  well  from 
earthquakes  in  the  same  regions  as  the  explosions. 
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I.  INTRODUCTION 

In  this  report  we  describe  a study  of  the  variability  of  the  shape  of 
the  amplitude  spectrum  In  the  main  signal  bandwidth  of  short-period  P-wave 
teleseisms  as  seen  on  LASA  recordings.  The  liiq>etus  for  such  a study  comes 
from  the  numerous  published  suggestions,  based  on  both  theory  and  observa- 
tion, that  ratios  of  certain  parts  of  this  P— wave  signal  band  may  effectively 
serve  as  a means  of  distinguishing  source  type,  l.e.,  underground  e3q>loslon 
or  natural  earthquake.  Before  explaining  the  purpose  and  procedure  of  our 
study  further.  It  seems  appropriate  to  review  past  work  on  this  spectral  dis- 
criminant. 


Early  evidence  of  difference  In  P-wave  spectra  between  earthquakes  and  | 

underground  explosions  was  presented  by  Willis  (1963)  and  Frantti  (1963),  * 

I 

but  these  contributions  represented  only  a few  seismic  events  and  suggested 
no  quantitative  discriminant  parameter.  In  reports  by  Briscoe  (1966)  and 
Briscoe  and  Walsh  (1967),  a ratio  of  energy  In  two  spectral  bands  was  sug- 

I 

gested  as  such  a parameter  and  was  evaluated  for  short-period  P-waves  record- 
ed at  LASA  from  several  Eurasian  events.  They  used  the  spectral  bands  of  ^ 

i 

0.3-0. 7 Hz  and  1. 5-1.9  Hz.  By  analyzing  a large  population  of  Eurasian  j 

events  which  Included  many  explosions,  Kelly  (1968)  and  Lacoss  (1969)  brought  j 

the  short-period  spectral  ratio  Into  prominence  as  a discriminant  which  cor-  j 
rectly  classified  roughly  90%  of  their  event  sample.  Based  on  calculations 
of  S/N  ratios  for  each  event,  Lacoss  (using  the  spectral  bands  0.45-0.95  Hz 


1 


Willis,  D.E.  (1963),  Comparison  of  seismic  waves  generated  by  different  types 
of  source.  Bull.  Seism.  Soc.  Am. , v.  53,  p.  965-978. 

Frantti,  G.E.  (1963),  Energy  spetra  for  underground  explosions  and  earth- 
quakes, Bull.  Seism.  Soc.  Am. . v.  53,  p.  997-1006. 

Briscoe,  H.W.  (1966),  Ratios  of  spectral  densities.  In  Semiannual  Technical 
Summary — Seismic  Discrimination,  30  June  1967,  Lincoln  Laboratory, 
Lexington,  Massachusetts. 


Briscoe,  H.W.  and  R.  walsh  (1967),  Ratios  of  spectral  densities.  In  Semian- 
nual Technical  Stannary — Seismic  Discrimination,  30  June  1967,  Lincoln 
Laboratory,  Lexington,  Massachusetts. 

Kelly,  E.J.  (1968),  A study  of  two  short-period  discriminants.  Technical  Mote 
1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 
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and  1.45-1.95  Hz)  detennlned  that  the  lASA  threshold  for  discrimination  by 


spectral  ratio  was  at  approxlmatly  ]ii|^*4.8  for  Eurasian  events.  This  magni- 
tude Is  often  quoted  as  the  conservative  limit  of  positive  Identification 
by  the  discriminant. 

The  P-wave  spectral  discriminant  has  continued  to  be  studied  In  various 
forms.  Perhaps  the  slnqplest  Is  a measurement  of  the  dominant  or  Initial 
period  of  the  P-wave  signal  In  the  time  domain  (Noponen,  1975;  Lacoss,  1969). 
A form  eiiq>loyed  by  Anglin  (1971  and  Welchert  (1971)  Is  the  third  moment  of 
frequency  (IMF)  of  the  amplitude  spectrum  between  0.32  and  5.0  Hz,  which 
rather  successfully  separated  a large  sample  of  Eurasian  explosions  and 
earthquakes.  A recent  Innovation  Is  the  plotting  of  autoregressive  parame- 
ters, which  are  fit  to  the  time  series  of  seismic  P waves,  by  Tjosthelm 
(1975).  Again,  this  technique  resulted  In  fair  separation  of  a large  sample 
of  Eurasian  events.  Shumway  and  Blandford  (1974)  applied  classical  linear 
discrimination  to  the  raw  event  spectra,  using  a learning  and  test  populat- 
ion of  Eurasian  events.  They  found  good  results  and  suggested  that  those 
explosions  idilch  failed  the  discriminant  were  cratering  events.  These 
workers  also  suggested  that  the  reason  for  the  success  of  the  discriminant, 
when  applied  to  Eurasian  explosions,  teas  the  cancellation  of  low-frequency 
P by  pP. 


Lacoss,  R.T.  (1969),  A large-population  LASA  discrimination  experiment. 
Technical  Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 

Noponen,  I.  (1975),  Compressions!  wave  power  spectrum  from  seismic 
sources — Final  scientific  report.  Institute  of  Seismology,  U.  of 
Helsinki,  Helsinki,  Finland. 

Lacoss,  R.T.  (1969),  A large-population  lASA  discrimination  experiment. 
Technical  Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 

Anglin,  F.M.  (1971),  Discrimination  of  earthquakes  and  explosions  using 
short-period  seismic  array  data.  Nature,  v.  233,  p.  51-52. 

Welchart,  D.H.  (1971),  Short-period  spectral  discriminant  for  earthquake- 
explosion  differentiation,  Z,  Geophys . « v.  37,  p.  147-152. 

Tjosthelm,  D.  (1975),  Autoregressive  representation  of  Seismic  P-wave  sig- 
nals with  an  application  to  the  problem  of  short-period  discriminants 
Geophys.  v.  43,  p.  269-291. 
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to  the  very  narrow  spectral  range  usually  Included  In  a spectral  ratio  defini- 
tion. In  practice.  If  only  a small  band  Is  used,  the  true  character  of  real 
seismic  sources  can  overwhelm  the  characteristics  predicted  from  ideal  sources. 
Examples  of  such  are:  1)  many  earthquakes  are  undoubtedly  accompanied  by  sub- 
earthquakes (Blandford,  1975)  which  distort  the  spectrum,  2)  the  source-time 
function  of  earthquakes  may  be  very  irregular,  3)  the  explosion  spectra  will 
have  varying  shape  due  to  different  P-pP  intervals,  4)  pP  may  not  be  perfect- 
ly reflected  due  to  spallation  over  explosions,  and  5)  changes  in  detonation 
medium  and  decoupling  may  greatly  affect  the  explolon  comer  frequencies. 

The  brief  mention  of  these  theoretical  source  aspects  indicates  that 
some  scatter  in  spectral  ratios  should  be  expected  due  to  departures  from 
ideal  source  conditions.  But  the  chief  factor  responsible  for  scatter  in 
spectral  ratios  may  be  the  variation  of  Q throughout  the  earth.  That  this 
variations  exists  radially  from  the  center  of  the  earth  presents  little  pro- 
blems in  treating  spectral  ratio  of  P-waves;  but  the  well-documented  lateral 
variations  (Barazangl  et  al.,  1975;  Molnar  and  Oliver,  1969;  Burton,  1972) 
must  affect  spectral  ratios  considerably  in  a more  random  fashion. 


I 

I 

I 
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In  this  report  we  will  document  the  variability  of  the  P-wave  spectral 
ratio  discriminant  by  using  a large  set  of  P-waves  recorded  at  the  Montana 
LASA.  Such  documentation  will  enable  us  to  grasp  the  real  scatter  in  this 
discrimination  parameter,  to  assess  whether  it  is  more  stable  over  small 
source  regions,  and  to  make  some  inferences  about  Q variations  in  the  mantle. 
Our  purpose  is  to  provide  a backdrop  of  data  against  which  the  probable  ef- 
fectiveness of  this  discriminant  used  on  a global  scale  could  be  Judged. 


i 


Blandford,  R.R.  (1975),  A source  theory  for  complex  earthquakes.  Bull.  Seism.  | 

Soc.  Am. . V.  65,  p.  1385-1406.  I 

Barazangl,  M. , W.  Penlngton,  and  B.  Isacks  (1975),  Global  study  of  seismic  1 

wave  attenuation  in  the  upper  mantle  behind  island  arcs  using  pP  waves, 
jJ.  Geophys . Res. , v.  80,  p.  1079-1092.  . j 

i 

Molnar,  P.,  and  J.  Oliver  (1969).  Lateral  variations  of  attenuation  in  the  I 

upper  mantle  and  discontinuities  in  the  lithosphere,  Geophys . Res. , ^ J 

V.  74,  p.  2648-2682.  ; 


-10- 


t J 


i 


DATA 

The  P-wave  signals  In  this  study  were  those  reported  in  the  LASA  daily 
bulletin  prepared  by  the  SDAC.  A cutoff  of  roughly  95”  distance  is  used  for 
the  LASA  bulletin  report  on  epicenters.  The  signals  processed  were  recorded 
In  the  time  period  of  January  through  August  of  1972  and  number  about  <000  or 
on  the  order  of  15  per  day.  The  signals  were  recovered  from  LASA  Event  Pro- 
cessor (EP)  tapes  by  association  with  the  LASA  bulletin.  We  recovered  only 
the  AO  subarray  beam  on  the  signal  to  avoid  signal  loss  which  occurs  at  high 
frequencies  In  the  full  array  beam  since  the  signal  correlation  length  decre- 
ases with  Increasing  frequency.  These  signals  were  already  screened  by  the 
LASA  Detection  Processor  (DP)  i^lch  allows  only  detections  (on  the  full  array 
beam)  with  a signal- to-nolse  ratio  > 5(14  dB)  to  be  recorded  on  the  EP  tapes. 
Thus  our  data  base  consisted  of  signals  which  already  had  fair  S/N  ratio,  as- 
suming that  AO  was  not  systematically  low  In  amplitude  relative  to  the  full 
array  beam.  Further  S/N  tests,  to  be  described  later  In  this  report,  were 
made  on  the  AO  subarray  signals  In  our  analysis  of  their  spectral  content. 

The  sampling  rate  of  the  data  was  10  samples/second. 

Hypocenters  used  In  this  report  are  those  reported  In  the  LASA  bulletin. 
Considerable  error  In  the  epicenters  exists  due  to  the  poor  location  capabil- 
ity of  a single  array;  this  error  has  been  estimated  theoretically  to  have  a 
standard  deviation  of  2.5*  by  Shllen  and  Toksoz,  1973.  Actual  comparisons 
with  NEIS  epicenters  (Shllen  and  Toksoz,  1973;  Ahner,  1973;  and  Woolson,  1976) 
show  that  the  standard  deviation  is  roughly  2*-3*.  Depths  reported  in  the 
LASA  bulletin  are  based  on  analysts'  decisions  on  the  appearance  of  pP  In  the 
LASA  record;  since  these  picks  are  uncorroborated  by  those  from  other  stations 
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•Shllen,  S.  and  M.N.  Toksoz  (1973),  Automatic  event  detection  and  location 
capabilities  of  large  aperture  seismic  arrays.  Bull.  Seism.  Soc.  Am. . 
V.  63,  p.  1275-1288. 

A 

Ahner,  R.O.  (1973),  A comparison  of  the  location  refinement  techniques  In 
the  SDAC  LASA  event  processor.  Report  No.  TR-73-5,  Seismic  Data  An- 
alysis Center,  Teledyne  Geotech,  Alexandria,  Virginia. 

Woolson,  J.R.  (1976),  unpublished  memorandum.  Teledyne  Geotech,  Alexand- 
ria, Virginia. 
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or  by  wiltlstatlon  hypocenter  location,  they  are  not  highly  reliable.  No  coo- 
paratlve  studies  were  available  with  which  to  quantitatively  Judge  the  accur- 
acy of  the  reported  depths  of  focus  In  the  LASA  bulletin  though.  Furthemore, 
■any  pP  are  undoubtedly  missed,  and  a normal  (33  km)  depth  Is  arbitrarily  as- 
signed by  the  analyst. 


Since  both  the  epicenter  and  depth  estimates  from  the  LASA  bulletin  are 
used  in  presenting  various  results  later  In  this  report,  we  feel  that  the  er- 
rors associated  with  them  will  have  some  adverse  l]iq>act  on  those  results. 

This  Is  especially  true  for  depth  errors  because  deep  earthquakes  may  have 
systematically  different  spectra  than  shallow  earthquakes.  An  alternative 
would  have  been  to  collate  our  LASA  signals  with  NETS  hypocenters,  but  the 
effort  Involved  In  that  procedure  did  not  seem  justified  In  relation  to  what 
we  felt  would  be  only  a small  gain  In  the  clarity  and  Importance  of  the  re- 
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METHODS  OF  CALCULATION 


Signals  taken  from  the  EP  tapes  are  shaped  by  the  LASA  Instrument  re- 
sponse which  Is  highly  peaked  at  4 Hz,  and  we  did  not  remove  this  response  In 
out  processing  of  the  data  since  we  are  concerned  with  relative  spectral  val- 
ues, Figure  1 Illustrates  the  output  of  our  processing  for  typical  P-waves 
at  LASA.  For  each  of  the  events,  the  four  traces  represent  outputs  of  filters 
with  bandpasses,  respectively,  of  0. 4-3.0,  0.4-0. 8,  1.4-1. 8,  and  1. 8-4.0  Hz. 
These  filters  are  applied  In  the  frequency  domain  by  simulating  the  responses 
of  analog  filters  with  these  comer  frequencies  and  with  rolloffs  of  24  db/oct 
on  either  side.  The  amplitude  of  the  traces  Is  equalized  for  plotting,  with 
the  maximum  amplitude  In  counts  given  at  the  left.  Because  of  the  peaked  In- 
striment  response,  the  filter  outputs  do  not  necessarily  have  their  predomin- 
ant frequencies  near  the  center  of  the  filter  bandpass. 

Within  a four— second  window  beginning  at  the  maximum  aiiq>lltude  In  the 
0. 4-3.0  Hz  filtered  P-wave  as  shown,  the  maximum  amplitude  (counts)  In  each  of 
the  four  filter  outputs  for  an  event  was  determined.  These  time-domain  ampli- 
tudes are  the  estimates  for  spectral  amplitude  within  each  passband.  Previous 
work  on  spectral  ratios  has  employed  frequency-domain  estimates  exclusively. 
But  since  these  estimates  were  usually  based  on  sample  lengths  of  10  or  even 
as  much  as  20  seconds  of  data  following  P onset,  there  was  not  assurance  that 
the  spectral  content  of  the  direct  P-wave  alone  was  accurately  represented  In 
the  conq}uted  spectra. 

The  S/N  ratio  for  each  passband  was  computed  by  forming  the  ratio  of 
maximum  amplitude  In  the  four-second  signal  window  to  the  maximum  amplitude 
within  the  noise  shown  preceding  the  P-wave.  These  S/N  ratios  and  the  "spec- 
tral" amplitudes  were  then  stored  along  with  the  LASA  hypocenter  data.  For 
the  results  presented  In  the  remainder  of  this  report,  any  bandpass  output 
from  the  AO  subarray  with  S/N  < 2 was  not  used.  This  requirement  more  than 
halved  the  number  of  available  signals. 


RESULTS 


Results  for  the  LL  Ratio 

In  this  report  we  define  the  "LL"  spectral  ratio  as  the  ratio  of  maxl- 
Bum  amplitude  In  the  0.4-0. 8 Hz  bandpass  output  to  that  In  the  1.4-1. 8 Hz 
output.  This  Is  a variant  of  the  original  Briscoe  (1966)  spectral  ratio,  and 
It  Is  essentially  the  Inverse  of  that  presented  In  many  reports  by  the  Appli- 
ed Seismology  Group  at  Lincoln  Laboratory. 


Since  our  purpose  Is  to  examine  the  variation  of  this  spectral  ratio 

over  different  regions  of  the  earth.  It  Is  desirable  to  Investigate  whether 

any  systematic  effects  are  present  In  the  computed  ratios.  One  effect  that 

follows  from  source  theory  Is  that  the  spectra  within  our  frequency  band 

should  be  a function  of  source  spatial  extent,  or  magnitude.  If  the  stress 

drop  Is  fixed.  Most  source  theories  predict  a teleselsmlc  pulse  with  a flat 

"2  —3 

low-frequency  spectrum  falling  off  as  f or  f above  a comer  frequency 

_2 

which  decreases  trlth  Increasing  magnitude.  For  a f model,  the  difference 

In  the  LL  spectral  ratio  between  very  low  and  very  high  magnitude  events 

should  be  the  difference  between  amplitude  ratios  made  at  the  center  freque- 

_2 

ncles  of  0.6  and  1.6  Hz  for  a flat  (low  m^^)  and  a f (high  nij^)  spectral  shape 
respectively.  This  amounts  to  a factor  of  about  seven.  Thus  the  LL  ratio  as 
defined  above  should  be  nearly  an  order  of  magnitude  greater  for  large  events 
compared  to  small  events.  In  figure  2 are  presented  all  LL  ratios  from  our 
data  base  for  which  the  S/N  ratio  exceeded  2.0  In  both  bandpass  outputs.  Also 
only  shallow  events  (<60  km,  as  reported  on  LASA  bulletin)  were  used.  This 
figure  shows  no  apparent  trend  for  the  LL  ratio  with  Increasing  m^^,  except  be- 
yond 5.5  perhaps;  apparently  In  disagreement  with  the  above  argument.  One  rea- 
son for  the  real  data  not  to  follow  theoretical  trends  may  be  that  the  small 
sample  of  events  above  mj^~5  Is  biased.  Another  reason  may  be  that  the  Instru- 
ment response,  peaked  at  3-4  Hz,  distorts  the  time-domain  estimate  of  spectral 
amplitude  after  application  of  narrow-band  filters  by  leakage  of  higher  freq- 


Brlscoe,  H.  W.  (1966),  Ratios  of  spectral  densities.  In  Semiannual  Technical 
Summary — Seismic  Discrimination,  31  December  1966,  Lincoln  Laboratory 
Lexington,  Massachusets . 
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usncy  coaponents  Into  lower  bandpass  outputs*  A check  of  several  events  such 
as  those  in  Figure  1 revealed  that  the  dominant  periods  fdLthln  the  signal  %rLn- 
dow  for  the  second  and  third  traces  were  always  close  to  the  center  frequen- 
cies of  0.6  and  1.6  Hz  though,  and  so  this  latter  reason  seems  precluded.  A 
remaining  possibility  Is  that  the  saaple  of  ratios  from  lower  magnitude  events 
are  heavily  weighted  with  short  eplcentral  distances  and  that  a bias  arises 
d\ie  to  different  attenuation  effects  between  near  and  far  events.  This  pos- 
sibility can  be  checked  by  plotting  the  data  of  Figure  2 versus  eplcentral  dis- 
tance as  In  Figure  3.  Here  %ie  have  restricted  the  data  base  to  m|^  > 4.6  In 

order  to  Insure  essentially  complete  detection  of  events  out  to  A > 95*  for 

the  AO  subarray  beam.  This  mj^  cutoff  Is  based  on  a full  array  90%  threshold 
of  % 3.9  (Dean,  1971)  iq>  to  A ■ 80*  and  a beaming  loss  of  vif  In  the  S/N 

ratio  where  M ■ 25  In  going  from  the  full  LASA  array  composed  of  25  subarrays 

to  the  AO  subarray.  Although  considerable  scatter  exists  In  the  data  of  Fig- 
ure 3,  there  Is  no  Indication  of  a trend  with  distance;  and  Inference  of  a near- 
ly constant  attenuation  factor  Is  not  unreasonable.  This  Implies  constant  t* 

In  the  equation 

A^(f)  - Ag(f)  e “’'^*'*.G(A) 

where  subscripts  o^  and  ^ mean  observed  and  source  spectrum  and  G(A)  Is  the 
geometrical  spreading  factor  which  Is  assumed  to  be  Independent  of  frequency. 

The  parameter  t*  Is  given  by 

t*  - / 

■'s  oQ 

frtiere  a Is  the  compresslonal-wave  velocity,  Q Is  the  quality  factor  for  com- 
presslonal  waves,  and  £ Is  the  path  length. 

Thus,  we  eliminate  the  possibility  of  significant  distance  bias  In  the 
results  of  Figure  2.  The  lack  of  empirical  trend  with  m^  for  the  LL  ratio  Is 
then  unexplained,  and  we  can  only  conclude  that  some  unapparent  bias  In  the 
data  exists  or  that  there  Is  really  a negligible  average  change  with  iBj^.  The 
negligible  change  could  more  easily  be  accommodated  by  the  f ^ spectral  slope 
at  Intermediate  frequencies  suggested  by  Brune  (1970)  In  the  case  of  partial 

Dean,  H.C.  (1971),  Detection  threshold  of  the  LASA/SAAC  system.  Report  No.  3, 

Seismic  Array  Analysis  Center,  Teledyne  Geotech,  Alexandria,  Virginia. 
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EPICENTRAL  DISTANCE,  DEGREES 

LL  spectral  ratio  versus  eplcentral  distance  from  LASA  for 
events  with  m,  >4.6  and  DOF  < 60  km. 


stress  drop,  which  is  probably  the  rule  rather  than  the  exception  for  earth- 
quakes. Another  possibility  Is  that  numerous  subearthquakes  (Blandford, 

1975)  make  the  high  frequency  spectrum  have  a smaller  slope.  Still  another 
possibility  Is  that  for  each  fixed  m^  the  partial  stress  drops  for  all  earth- 
quakes at  that  are  widely  varying.  Thus  at  each  mj^  there  would  be  a wide 

range  of  comer  frequencies,  obscuring  a trend  toward  more  low  frequencies 

* 

at  high  magnitudes. 


Together  Figures  2 and  3 suggest  that  no  systematic  effects  are  contain- 
ed In  our  computed  spectral  ratios.  If  these  ratios  are  considered  on  a re-  ^ 
glonal  basis,  then  any  variation  from  region  to  region  should  be  reflective 
of  real  difference  either  In  source  spectral  character  or  In  the  parameter 
t*  along  the  particular  paths  to  LASA.  In  Figure  4 are  shown  the  values  > 

10«log^Q(LL)  where  log^QLL  Is  the  average  of  the  logarithm  of  the  spectral  i 

ratios  for  events  within  given  geographic  regions,  as  configured  by  Fllnn  et 
al.  (1974).  Only  regions  within  roughly  100**  of  LASA  are  represented  due  to 
the  LASA  bulletin  cutoff,  and  one  would  not  want  to  consider  more  distant 
regions  an}way  since  effects  of  the  core  are  present  on  P-waves  from  eplcen-  I 
ters  beyond  roughly  95**.  Many  regions  produced  no  data  for  this  plot,  and 
also  regions  with  only  one  signal  were  Ignored.  A plot  of  the  166  geograph- 
ic region  numbers  represented  In  Figure  4 Is  shown  In  Figure  5 for  refer- 
ence; the  region  numbers  are  plotted  near  the  centroid  of  each  region.  We 
note  that  for  region  329  In  Central  Asia  the  value  plotted  In  Figure  4 Is 
based  on  four  presumed  underground  explosions  and  cannot  be  coi^>ared  to 
earthquake  results.  Table  I lists  the  data  used  In  plotting  Figure  4 and 
also  the  sample  size  NEQ  and  standard  deviation  SD  of  the  Individual  data 
for  each  region.  The  approximate  confidence  limits  for  the  mean  log^^CLL) 

In  each  region  Is  given  by  +(2*SD)/  »^EQ,  and  It  Is  apparent  that  much  of 
the  variation  seen  In  Figure  4 Is  well  established  by  statistical  measure. 
Examination  of  a map  of  SD  values  similar  to  Figure  4 but  not  shown  here 
revealed  no  systematic  variation  clearly  associated  with  tectonic  character. 


Blandford,  R.R.  (1975),  A source  theory  for  complex  earthquakes.  Bull.  Slesm. 
Soc.  Am. . V.  65,  p.  1385-1406. 

Fllnn,  E.A.,  E.R.  Engdahl,  and  A.R.  Hill  (1974),  Seismic  and  geographical  re- 
gionalization, Bull.  Seism.  Soc.  Am. . v.  64,  p.  771-992. 
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TABLE  I (Continued) 

Regionalized  Data  for  the  LL  Spectral  Ratio 


log^^Q  (LL  spectral  ratio) 


CRN 

MEAN 

ST. DEV. 

NEC 

210 

.167 

.198 

7 

211 

- .013 

.254 

73 

212 

-091 

.223 

10 

213 

.097 

.136 

11 

215 

.113 

.194 

9 

216 

.181 

.305 

7 

217 

.109 

.232 

9 

218 

- .107 

.213 

20 

219 

.022 

.201 

20 

220 

- .277 

.053 

3 

221 

- .107 

.197 

42 

222 

.011 

.209 

11 

223 

- .168 

.262 

4 

224 

- .013 

.208 

10 

226 

- .144 

.274 

2 

227 

.018 

.147 

5 

228 

- .139 

.328 

12 

229 

.051 

.205 

41 

232 

- .042 

.024 

3 

233 

- .127 

.197 

7 

234 

.473 

.300 

3 

235 

.297 

.365 

3 

236 

.126 

.454 

2 

238 

.207 

.037 

2 

245 

.331 

.181 

12 

321 

.019 

.209 

17 

323 

.349 

.308 

3 

329 

- .288 

.146 

4 

332 

.339 

.223 

7 

334 

.188 

.414 

2 

344 

.099 

.158 

4 

357 

- .288 

.156 

2 

362 

.309 

.056 

3 

363 

.006 

.098 

6 

364 

.299 

.261 

8 

365 

.136 

.230 

3 

366 

.195 

.226 

21 

368 

.438 

.166 

3 

370 

.120 

.181 

6 

371 

.262 

.062 

5 

375 

.078 

.099 

4 

381 

.246 

.342 

4 

24 


TABIiE  I (Continued) 

Regionalized  Data  for  the  LL  Spectral  Ratio 
lo^io  (IJ*  spectral  ratio) 


CRN 


MEAN 


ST.DF.V. 


-25- 


NEQ 


382 

.013 

.146 

2 

390 

.208 

.139 

3 

400 

.144 

.297 

4 

402 

.243 

.247 

47 

403 

.344 

.206 

24 

404 

.219 

.161 

6 

406 

.234 

.246 

32 

515 

.036 

.030 

2 

525 

.290 

.275 

8 

543 

.286 

.140 

2 

545 

.016 

.041 

3 

611 

.133 

.144 

3 

612 

.228 

.069 

4 

614 

.324 

.144 

4 

618 

.358 

.008 

2 

633 

.074 

.231 

7 

635 

.097 

.220 

4 

636 

.446 

.325 

5 

637 

.390 

.199 

12 

638 

.636 

.129 

2 

640 

.575 

.257 

6 

656 

.145 

.261 

5 

658 

.004 

.075 

2 

659 

- .145 

.357 

4 

661 

- .117 

.016 

2 

662 

.150 

.436 

2 

663 

- .176 

.198 

7 

666 

.496 

.230 

2 

670 

.106 

.169 

2 

677 

.174 

.150 

4 

684 

.293 

.218 

9 

685 

.262 

.264 

16 

686 

.252 

.023 

2 

692 

.320 

.319 

6 

693 

.354 

.163 

18 

694 

.372 

.278 

26 

695 

.572 

.255 

11 

696 

.260 

.311 

2 

697 

.348 

.210 

2 

715 

- .201 

.206 

7 

717 

- .175 

.000 

2 

724 

- .262 

.090 

3 

Results  for  the  HF  Ratio 

Before  interpreting  the  results  contained  in  Figure  4,  we  will 
another  spectral  ratio  defined  using  our  highest  bandpass  filter  froa  1.8- 
4.0  Hz.  The  ratio  of  the  maximum  aiiq)lltude  in  the  0.4-0. 8 Hz  output  to  that 
in  this  output  is  termed  here  the  "HF"  spectral  ratio.  Presumably,  if  Q ef- 
fects are  highly  variable,  this  ratio  should  have  even  greater  scatter  than 
the  LL  ratio.  Figure  6 is  a plot  of  this  data,  with  the  restriction  of  S/H 
ratio  > 2.0  for  both  traces  as  for  the  LL  ratio.  Comparison  of  Figure  6 %rlth 
Figure  2 indicates  no  significant  Increase  in  scatter,  and  so  effects  other 
than  Q on  the  spectra  must  be  large.  The  fitted  line  in  Figure  6 does  indi- 
cate a change  in  spectral  ratio  with  m^,  as  predicted  from  source  theory,  but 
not  nearly  as  large  as  expected. 

Figure  7 shows  the  data  of  Figure  6,  with  a cutoff  below  m^  « 4.6,  plot- 
ted versus  epicentral  distance;  and,  as  for  the  LL  ratio  in  Figure  3,  no 
systematic  trend  with  distance  is  apparent  for  the  HF  ratio. 

Figure  8 presents  the  regional  averages  for  the  logarithms  of  the  HF 
ratio,  and  Table  II  lists  these  data.  Only  134  regions  are  represented  here, - 
32  less  th£m  for  the  LL  ratios  in  Figure  4,  because  many  more  1. 8-4.0  Hz 
bandpass  outputs  were  screened  out  by  the  S/N  ratio  requirement  than  1.4-1. 8 
hz  outputs.  The  regional  pattern  shown  in  Figure  4 is  essentially  repeated 
by  Figure  8. 

Q Variations 

Neglecting  all  other  contributions  to  the  scatter  in  the  data,  let  us 
focus  on  the  possible  range  in  Q,  or  rather  t*,  if  this  parameter  alone  caused 
the  observed  variation  in  the  spectral  ratios.  (Later  we  will  show  why  t* 
alone  cannot  explain  the  scatter.)  By  using  equation  (1)  above,  the  differ- 
ence between  two  observed  values  of  LL  would  be 

(logLLj^-logLL2)  - iT(f^-fj^)  (tj^*-t2*)/loggl0  ^2) 

where  fjj  and  f^  can  be  taken  as  the  center  frequencies  of  the  two  passbands 
used  in  obtaining  LL,  namely,  0.6  and  1.6  Hz.  Using  these  frequency  values 
and  the  maximum  spread  in  log^^CLL)  of  roughly  1.6  from  Figure  2 for  the  left- 
hand  side  of  (2),  we  can  estimate  the  maximum  t*  spread  as  roughly  (tj^*-t2*) 
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TABLE  II 


Regionalized  Data  for  the  HF  Spectral  Ratio 
log^^Q  (HP  spectral  ratio) 


CRN 

MEAN 

ST. DEV. 

m. 

1 

.052 

.209 

7 

2 

.132 

.315 

6 

4 

.175 

.061 

7 

5 

- .044 

.133 

3 

6 

- .022 

.340 

9 

7 

- .031 

.280 

41 

9 

.015 

.258 

20 

10 

.155 

.158 

8 

12 

.047 

.379 

2 

13 

.184 

.401 

2 

15 

.127 

.192 

6 

16 

.015 

.497 

2 

17 

.275 

.252 

5 

19 

.626 

.195 

3 

20 

.251 

.337 

4 

21 

.254 

.635 

2 

46 

.005 

.152 

2 

49 

.090 

.180 

7 

51 

.503 

.001 

2 

53 

.111 

.179 

4 

60 

.264 

.292 

9 

61 

.276 

.266 

17 

62 

.234 

.228 

7 

63 

.118 

.368 

5 

65 

.085 

.052 

2 

66 

.263 

.023 

2 

69 

.266 

.304 

5 

70 

.200 

.245 

4 

71 

.349 

.340 

8 

74 

.082 

.096 

9 

75 

.276 

.329 

4 

76 

.206 

.271 

15 

78 

- .237 

.101 

2 

81 

.176 

.149 

2 

83 

.207 

.356 

9 

90 

- .018 

.389 

2 

92 

.164 

.422 

7 

94 

.149 

.340 

5 

95 

- .028 

.075 

3 

99 

- .303 

.224 

18 

103 

- .096 

.267 

13 

109 

.251 

.423 

4 
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TABLE  II  (Continued) 

Regionalized  Data  for  the  HF  Spectral  Ratio 
log^^Q  (HF  spectral  ratio) 


CRN 

MEAN 

ST.DE 

no 

.292 

.377 

111 

.054 

.167 

115 

.136 

.240 

116 

.063 

.338 

118 

- .030 

.156 

120 

.120 

.474 

122 

.163 

.292 

123 

.298 

.213 

124 

- .121 

.127 

125 

- .190 

.235 

127 

.124 

.509 

128 

- .072 

.048 

129 

- .055 

,170 

130 

.127 

.484 

134 

.475 

.212 

135 

.252 

.230 

136 

.376 

.336 

137 

.014 

.269 

146 

.108 

.102 

169 

.442 

.031 

170 

.597 

.002 

171 

.129 

.136 

173 

.232 

.282 

174 

- .063 

.496 

175 

.046 

.008 

181 

.282 

.263 

183 

.172 

.318 

210 

.385 

.166 

211 

.143 

.228 

212 

.243 

.270 

213 

.232 

.184 

215 

.380 

.304 

216 

.340 

.235 

217 

.286 

.268 

218 

.044 

.226 

219 

.226 

.265 

220 

- .080 

.112 

221 

.037 

.218 

222 

.184 

.265 

224 

.028 

.223 

226 

- .031 

.397 

227 

.206 

.155 

f ( 


1 


TABLE  II  (Continued) 

Regionalized  Data  for  the  HF  Spectral  Ratio 
log^^Q  (HF  spectral  ratio) 


I 


CRN 

MEAN 

ST. DEV. 

m. 

228 

.120 

.263 

10 

229 

.204 

.215 

30 

230 

.000 

.180 

2 

232 

.191 

.128 

4 

233 

.253 

.313 

8 

235 

.158 

.187 

2 

245 

.475 

.268 

9 

321 

.107 

.145 

10 

323 

.516 

.240 

2 

329 

- .241 

.221 

4 

332 

.347 

.343 

4 

344 

.389 

.091 

2 

357 

- .181 

.017 

2 

363 

.378 

.136 

5 

364 

.358 

.389 

5 

365 

.149 

.067 

2 

366 

.144 

.209 

10 

368 

.502 

.111 

3 

370 

.384 

.208 

5 

381 

.386 

.337 

6 

400 

.516 

.333 

2 

402 

.422 

.273 

32 

403 

.379 

.241 

13 

404 

.563 

.155 

5 

406 

.345 

.250 

14 

525 

.354 

.136 

5 

611 

.366 

.085 

3 

612 

.557 

.137 

2 

633 

.571 

.255 

4 

635 

.506 

.332 

3 

636 

- .221 

.467 

3 

637 

.320 

.502 

2 

640 

.185 

.245 

2 

656 

.232 

.423 

5 

657 

.190 

.558 

2 

658 

.095 

.051 

2 

659 

.354 

.121 

2 

661 

.024 

.000 

2 

662 

.386 

.337 

3 

663 

- .045 

.381 

5 

677 

.441 

.261 

3 

i 
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TABLE  II  (Continued) 

4P 

1 

Regionalized  Data  for  the 

HP  Spectral 

Ratio  ; 

(HP  spectral 

ratio) 

GRN 

MEAN 

ST. DEV. 

NE2 

682 

.214 

.007 

2 1 

683 

.356 

.099 

2 

684 

.353 

.162 

5 

685 

.393 

.176 

9 

692 

.497 

.185 

3 

693 

.636 

.157 

5 

694 

.186 

.324 

12 

696 

.403 

.034 

2 

715 

- .106 

.000 

2 
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^ 1.2.  A similar  calculation  for  the  HF  ratio  using  a spread  of  1.6  from 
Figure  6 and  fjj  ••  2.9  hz  gives  roughly  (tj^*-t2*)  ^ .5.  However,  the  data  in 
Figure  6 do  not  Include  many  events  which  appeared  In  Figure  2 (1022  versus 
1438)  because  the  higher  bandpass  output  did  not  pass  the  S/N  test.  The 
ratios  In  these  cases  would  generally  scatter  towards  the  top  of  the  measur- 
able HF  ratios  shown  In  Figure  5 and  probably  Increase  the  total  scatter  sig- 
nificantly If  there  had  been  S/N  > 2.0  on  the  1. 8-4.0  Hz  output.  Also  the  high 
frequency  band  probably  cannot  be  characterized  by  its  center  frequency  of 
2.9  hz  since  the  spectrum  changes  so  rapidly  In  the  band  1. 8-4.0  Hz.  Thus, 
we  regard  At*  £ 1.2  estimated  from  the  LL  ratios  as  more  reliable.  If  we 
assume  a minimum  possible  t*  of  0.2  for  a very  high  Q path  to  LASA  (Der, 

1976),  then  the  maximum  t*  Is  1.4,  giving  a range  which  is  not  unacceptable. 
However,  we  note  that  much  of  the  scatter  In  Figure  2 must  be  source  related 
(pP  Interference,  comer- frequency  variation,  time-function  variation,  etc.); 
and  so  the  Inferred  range  of  t*  Is  surely  less  than  1.2. 

Scatter  and  Regional  Variation  In  Spectral  Ratios 

In  relation  to  discrimination,  the  outstanding  feature  of  the  data  pre- 
sented here  Is  the  large  scatter  In  the  spectral  ratios  as  best  seen  In  Fig- 
ures 2 and  6.  Roughly  the  same  scatter  can  be  found  In  the  LL  ratio  for  156 
Eurasian  earthquakes  analyzed  by  Lacoss  (1969).  This  scatter  ranges  up  to 
1.6  logarithmic  units  and  is  more  than  half  as  large  as  that  in  obser- 

ved at  LASA  at  a given  magnitude  for  earthquakes  from  several  regions  (Capon 
et  al.,  1969).  This  occurs  for  the  LL  ratio  even  though  we  have  Included 
only  a very  small  spectral  bandwidth  In  the  estimate  of  this  ratio,  0.4  to 
1.8  Hz.  The  scatter  in  LL  and  HF  ratios  can  be  explained  by  acceptable  var- 
iations In  t*,  or  Q,  over  the  various  paths  to  LASA  although  many  source- 
related  effects  must  contribute  also. 

It  Is  reasonable  to  presume  that  these  spectral  variations  are  coherent 


Lacoss,  R.T.  (1969),  A large-population  LASA  discrimination  experiment. 

Technical  Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 

Capon,  J.,  R.J.  Greenfield,  and  R.T.  Lacoss  (1960),  Long-period  signal 

processing  results  for  the  Large  Aperature  Seismic  Array,  Geophysics « 
V.  34,  p.  305-329. 


over  some  range  of  azimuth  and  distance  from  LASA,  l.e.,  are  regional  In  cha- 
racter. Figures  4 and  8 can  be  Interpreted  as  a verification  of  this  pre- 
sumption, especially  since  examination  of  these  maps  reveals  a similarity  be- 
tween the  location  of  high  and  low  attenuating  source  regions  defined  by  the 
spectral  ratios  and  the  location  of  such  regions  from  other  types  of  studies 
(P  and  pP  waves  by  Barazangl  et  al.,  1975;  S waves  by  Molnar  and  Oliver, 

1969) . Oceanlc-rldge  source  regions  have  consistently  high  ratios  In  our 
figures.  This  suggests,  but  does  not  demand,  high  attenuation,  as  also  In- 
ferred by  Molnar  and  Oliver.  Due  to  the  large  errors  In  the  LASA  epicenters. 
It  Is  not  possible  to  clearly  distinguish  the  adjacent  high  and  low  attenua- 
tion zones  which  were  well  defined  by  Barazangl  et  al.  In  regions  of  plate 
subductlon.  Also,  most  of  the  events  east  of  Chile,  and  In  other  regions 
characterized  by  deep  events,  must  be  deeper  than  60  km;  but  pP  was  not  de- 
tected by  the  LASA  analyst,  which  In  Itself  suggests  higher  attenuation  In 
the  zone  above  the  earthquakes.  Since  our  method  of  analysis  was  not  desig- 
ned to  study  attenuation  directly,  only  crude  Inferences  can  be  made  from 
our  spectral  ratios. 

Because  much  of  the  overall  scatter  In  the  LL  or  HF  spectral  ratio  Is 
due  to  regional  variations,  this  scatter  could  be  reduced  by  regional  cor- 
rections as  has  often  been  suggested  for  discrimination.  Yet,  Tables 

1 and  11  reveal  that  considerable  scatter  still  exists  for  values  within  a 
given  region.  A typical  standard  deviation  Is  perhaps  0.20  for  the  loga- 
rithms of  a region's  ratios  while  that  for  all  the  ratios  In  both  Figures  2 
and  6 Is  0.30.  While  we  must  recall  that  the  errors  In  location  by  LASA 
will  cause  some  Intermixing  of  values  from  adjacent  regions,  we  fully  expect 
that  most  of  this  scatter  would  remain  with  highly  accurate  locations  since 
It  Is  due  to  variations  In  source  character  and  depth  of  focus  and  to  later- 
al Inhomogeneltles  within  a given  region.  Moreover,  If  the  parameters  M^ 

Barazangl,  M. , W.  Penlngton,  and  B.  Isacks  (1975),  Global  study  of  seismic 

wave  attenuation  In  the  upper  mantle  behind  Island  arcs  using  pP  waves. 
Geophysics.  Res. . v.  80,  p.  1079-1092. 

Molnar,  P.  and  J.  Oliver  (1969),  Lateral  variation  of  attenuation  In  the 

iq>per  mantle  and  discontinuities  In  the  lithosphere,  Geophys.  Res.. 

V.  74,  ?.  2648-2682. 
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and  uij^  are  any  guide,  subdividing  into  very  small  source  regions  may  not  gre- 
atly reduce  the  scatter  within  each  of  these  small  regions;  e.g.,  the  study 
of  LASA  in  three  Kamchatka  subregions  by  Blandford  and  Clark  (1975) 

showed  that  at  least  1.5  orders  of  magnitude  variation  was  present  in  at  a 
given  mjj  for  shallow  earthquakes  all  within  an  area  defined  by  a circle  with 
roughly  1**  diameter. 

Explosion  Ratios 

In  Figures  2 and  6 we  have  identified  those  spectral  ratios  for  seven 
presumed  imderground  nuclear  explosions  in  the  USSR.  The  epicenters  were 
identified  as  possible  explosions  by  one  or  more  of  these  criteria;  1)  as- 
signed zero  depth  on  the  NEIS  list,  2)  their  appearance  in  the  Uppsala,  Swe- 
den, reports  on  Soviet  explosions,  and  3)  their  low  values  in  von  Seggem 
(1976).  In  general  these  explosion  points,  either  the  LL  or  HF  ratios,  lie 
among  the  very  lowest  of  computed  values  for  all  the  earthquakes.  We  note 
again  that  our  data  is  probably  contaminated  with  a large  number  of  deep  ev- 
ents, but  not  reported  as  such  on  the  LASA  bulletin;  these  events  may  contri- 
bute heavily  to  the  lower  spectral  ratio  values  in  Figure  2 and  6.  We  must 
not  established  the  degree  of  attenuation  along  the  paths  from  the  Soviet 
test  sites  to  LASA  relative  to  other  paths.  In  our  data  base,  the  number  of 
natural  events  in  those  geographical  regions  which  Include  the  Soviet  explo- 
sion sites  in  Kazakh,  the  Urals,  and  above  the  Black  Sea  are  insufficient  to 
determine  t*  with  any  confidence.  From  spectral  analysis  of  explosions  though 
von  Seggern  and  Sobel  (1976)  have  determined  t*  for  Kazakh  to  LASA  to  be 
roughly  0.4,  a value  which  alone  cannot  account  for  the  appearsnce  of  the  ex- 
plosion points  so  low  within  the  global  sample  of  earthquakes. 


Blandford  R.R.  and  D.M.  Clark  (1975),  Variability  of  seismic  waveforms  at 
LASA  from  small  subregions  of  ECamchatka,  Report  No.  TR-75-12,  Seismic 
Data  Analysis  Center,  Teledyne  Geotech,  Alexandria,  Virginia. 

von  Seggem,  D.H.,  (1976).  Final  report  on  the  analysis  of  the  VLPE  network 
Report  No.  TR-76-1,  Seismic  Data  Analysis  Center,  Teledyne  Geotech, 
Alexandria,  Virginia. 

von  Seggem,  D.H.,  and  P.A.  Sobel  (1976).  Study  of  selected  Kamchatka 

earthquakes  in  a seismic  discrimination  context.  Report  No.  TR-76- , 

Seismic  Data  Analysis  Center,  Teledyne  Geotech,  Alexandria,  Virginia. 


In  Figure  9 we  present  the  explosion  LL  ratios  with  only  those  ratios 
from  earthquakes  in  a limited  area  composed  of  10  geographic  regions  aroimd 
the  USSR  explosion  sites.  These  regions  together  comprise  an  area  having  no 
known  rift  or  subductlon  zones,  except  for  the  Hindu  Kush  seismic  zone,  and 
which  should  therefore  be  rather  homogeneous  in  regard  to  Q properties.  The 
seven  explosion  points  do  not  separate  well  from  the  earthquake  population. 
Evidence  that  the  remaining  points  are  Indeed  earthquakes  and  not  explosions 
comes  from  their  values  (von  Seggem,  1976)  which  are  not  einomalous  and 

lie  well  above  those  of  the  seven  ejqploslons.  Also,  we  attempted  to  collate 
all  earthquakes  with  logLL  < 0 in  this  figure  with  the  NEIS  list.  The  col- 
lated ones  are  circled,  and  in  all  cases  the  NEIS  depth  was  < 70  km.  Thus 
low  earthquake  spectral  ratios  are  not  due  to  deep  Hindu  Kush  events. 

Corrections  to  m^  Based  on  Observed  Spectral  Ratios 


In  this  section,  an  attempt  will  be  made  to  utilize  the  average  spectral 
ratios  for  various  geographic  regions,  as  determined  In  this  report,  for  est- 
imating an  mj^  correction.  The  premise  that  the  spectral  ratios  are  an  Indi- 
cator of  relative  attenuation  has  already  been  shown  to  be  tenable  by  its 
correlation  with  known  geophysical  data.  There  are  several  factors  contri- 
buting to  the  spectral  ratio  parameter  and  to  m^  such  as  stress  drop,  comer 
frequency,  and  Incoherence  in  the  source  time  function;  nonetheless,  we  ex- 
pect that  our  spectral  ratio  results  can  be  helpful  in  reducing  the  scatter 
of  m^,  as  estimated  at  LASA,  for  a given  true  magnitude.  Since  the  true  mag- 
nitudes are  unavailable,  we  take  as  the  reference  magnitude.  We  assume 

without  evidence  that  variations  do  not  correlate  trlth  the  m^  variations. 

This  should  be  especially  correct  since  we  will  use  net%rork-determlned 
values  (von  Seggem,  1976)  for  which  LASA  was  not  a reporting  station. 

Mapping  the  spectral  ratio  Into  an  m^  correction  Is  done  by  first  using 
equation  (2)  to  obtain  the  t*  variation.  For  those  regions  which  had  enough 
events  to  calculate  an  average  logLL  as  given  in  Table  I,  we  see  that  the 
spread  Is  roughly  1.0  log  unit,  which  when  substituted  In  (2),  gives  a At* 


von  Seggem,  D.H.  (1976),  Final  report  on  the  analysis  of  the  VLPE  network. 
Report  No.  TR-76-1,  Seismic  Data  Analysis  Center,  Teledyne  Geotech, 
Alexandria,  Virginia. 
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spectral  ratios  for  events  In  geographic  regions  321,  329 
336,  344,  357,  362,  715,  717  and  724  (see  Figure  5). 


of  0.75.  Substituting  this  At*  value  in  equation  (1)  and  taking  logarithms 
results  in  an  Aiii|^  (assuming  £ - 1 Hz)  of  approximately  1.0  between  signals 
along  the  lowest  and  highest  attenuating  paths  to  LASA.  Thus,  we  propose 
this  simple  correction  to  LASA  m.  's: 


The  proposed  corrections  then  are  Identical  to  the  numbers  listed  In  Table  1 
We  will  use  only  those  regions  (108)  which  have  NEQ  > 4 In  this  table  how- 


The  LASA  in|^'s  were  collated  with  the  VLPE  values;  this  was  possible 
for  only  349  events  since  the  epicenter  lists  were  different  and  since  an 
was  not  reported  by  the  VLPE  network  for  many  events  within  the  95“  ring  sur- 
rounding LASA.  In  Figure  10a  are  plotted  separately  the  originally  reported 
LASA  mjj's  and  the  corrected  according  to  relation  (3)  versus  M^.  As 

clearly  represented  In  the  plots,  the  calculated  standard  deviation  of  the 
data  has  not  been  reduced.  The  standard  deviations  are  the  results  of  max- 
Imum-llkllhood  line  fitting  under  the  assumption  of  equal  variance  In  M 

8 

and  nij^.  Even  though  some  of  the  scatter  of  the  data  In  Figure  10a  Is  due  to 
M^,  that  part  due  to  mj^,  and  therefore  the  overall  scatter,  was  expected  to  be 
reduced.  This  result  suggests  that  the  spectral  ratios  have  little  actual 
correlation  with  attenuation  differences  Inferred  from  our  analysis  of  spec- 
tral ratios  versus  geographic  region.  To  affirm  this  result,  we  extracted 
only  those  events  from  geographic  regions  with  an  m^  correction  given  by  re- 
lation (3)  of  logLL  > 0.3  or  logLL  < -0.1.  This  left  only  110  events  from 
21  geographic  regions,  and  their  plot  are  shown  in  Figure  10b.  Again, 

even  though  we  have  taken  only  that  data  which  is  subject  to  large  correct- 
ions In  ntj^,  the  scatter  of  the  data  Is  not  significantly  reduced.  It  may  be 
that  have  over-corrected  the  m^  values  with  relation  (3);  to  test  this  Idea 
we  replotted  the  data  of  Figures  10a  and  10b  with  corrections  to  m^^  exactly 
one-half  those  already  used,  but  no  reduction  In  scatter  occurred. 


If  spectral  ratio  were  an  Indicator  of  attenuation  differences  and  thus 
variation  In  observed  m^  for  a given  source  spectrum,  then  we  would  fully  ex- 
pect reduction  In  m^  scatter.  In  fact,  we  have  not  demonstrated  any  l]iq>rove- 
ment;  and  this  Is  especially  unexpected  because  the  spectral  ratios  of  p-waves 
from  various  regions  seemed  to  correlate  well  with  direct  or  Indirect  geophys- 


ratio  for  the  geographic  reg: 


Figure  10b  Subset  of  data  from  Figure  10a  which  Includes  only  events  In 
those  geographic  regions  where  average  LL  ratio  was  > 0.3  or 


leal  evidence  on  lateral  Q variations  In  the  mantle.  An  alternative  explana- 
tion of  the  spectral  ratio  variations  would  be  the  source  spectrum  Itself. 

If,  due  to  systematic  change  In  stress  levels  or  other  physical  conditions 
affecting  earthquakes  mechanics,  there  were  regional  variations  In  source 
spectra  that  could  masquerade  as  Q effects,  then  these  variations  would  not 
necessarily  correlate  with  m^.  Another  alternative  e3q>lanatlon  would  be  that 
the  values  are  regionally  biased  such  that  corrections  to  m^^  only  serve  to 
move  the  point  farther  from  the  fitted  lines  In  Figures  10a  and  10b.  A third 
alternative  Is  that  the  Imprecision  of  the  spectral  ratio  as  a measure  of  at- 
tenuation or  the  errors  In  LASA  epicenters  have  so  affected  the  data  that 
reasonable  geophysical  hypotheses  cannot  be  confirmed. 


In  a study  of  wouldwlde  averages  from  NEIS  data,  Prozorov  and  Hud- 

son (1975)  have  shown  a pattern  In  this  parameter  which  Is  similar  to  our 
logLL  averages  shown  In  Figure  4,  and  they  Interpret  the  result  as  possibly 
due  to  attenuation.  Our  data  do  not  support  this  hypothesis,  and  evidence 
against  such  a conclusion  Is  the  plot  of  those  21  regions  having 

large  corrections  to  nij^,  as  used  In  Figure  10b.  Here,  VLPE  and  uncorrect- 
ed LASA  m^  are  used.  If  no  other  effects  were  present,  attenuation  would 
cause  the  parameter  to  Increase  with  logLL  as  Indicated  In  Figure  11; 

but  the  actual  data  scatter  greatly  and  do  not  support  such  a trend.  Figure. . 
11  clearly  shows  fdiy  our  attenq>ts  to  correct  m^  for  spectral  ratio  failed  to 
Improve  the  accuracy  of  the  m.  determination  at  LASA. 


1 

! 


Prozorov,  A.,  and  J.  A.  Hudson  (1974),  A study  of  the  magnitude  difference 
for  earthquakes,  Geophys.  ^. , v.  39,  p.  551-564. 
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Although  a simple  constant  stress-drop  source  theory  would  predict  a 
trend  with  magnitude  of  the  spectral  ratios  determined  in  this  report,  such 
a trend  is  barely  evident  and  is  much  less  than  would  be  predicted  theoret- 
ically. Either  some  unknown  biases  affect  the  data  or  more  complex  earth- 
quake spectral  models  are  required  to  satisfy  our  data.  Tentatively  we  sug- 
gest that  short-period  spectral  ratio,  at  least  formed  within  the  band  of 
0.4— 4.0  Hz,  can  be  studied  and  used  virtually  independent  of  magnitude  consid- 
erations. This  holds  also  for  distance  consideration  because  the  spectral 
ratios  were  found  to  be  constant  with  distance,  and  this  indicates  a rather 
constant  t*  versus  distance  relation. 

A total  logarithmic  scatter  in  our  LL  and  HF  spectral  ratios  was  found 
to  be  roughly  1.6  with  a standard  deviation  of  approximately  0.3.  These  val- 
ues are  consistent  with  previous  investigations  even  though  considerably  more 
source  regions  were  Included  in  this  study.  These  new  results  are  based  on 
computations  for  over  1400  earthquakes  teleselsmic  from  LASA,  a larger  data 
base  than  for  any  study  here-to-fore.  Such  scatter  is  less  than  that  obser- 
ved in  plots  for  global  events  but  represents  considerable  spectral 

variations  over  a small  bandwidth.  This  scatter  could  be  explained  by  effects 
on  the  spectra  of  reasonable  variations  in  the  t*  parameter,  from  roughly  0.2 
to  1.4;  but  we  favor  a smaller  range  in  t*  in  conjunction  with  several  source- 
related  effects  on  the  spectra  as  the  cause  of  the  observed  range  in  spectral 
ratios. 

Much  of  the  scatter  falls  into  definite  source-region  effects;  and  the 
resulting  pattern  of  regional  effects  bears  a similarity  to  known  upper  man- 
tle Q properties,  thus  implicating  Q variations  as  part  of  the  reason  for  the 
observed  scatter.  But  Intrareglonal  scatter  of  spectral  ratios  remains  high, 
with  a typical  value  of  0.2  for  the  standard  deviation,  and  may  be  caused  by 
variations  in  source  character.  This  implies  that  even  for  low  Q regions, 
the  probability  of  observing  an  "explosion-like"  spectral  ratio  is  not  negli- 
gible. For  high  Q regions,  the  probability  may  become  distressingly  high  in 
a test-ban  monitoring  scheme  which  used  spectral  ratio  as  a discriminant. 

This  extreme  variability  has  not  been  so  clearly  pointed  out  in  earlier  work 
by  other  authors. 
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In  this  report  we  attempted  for  the  first  time  to  reduce  scatter  In 
Mg-m^  plots  hy  jCorrectlng.  the  LASA  according  to  average  observed  spectral 
ratio  from  the, region  of  the  epicenter.  The  attempt  was  unsuccessful.  Our 
conclusion  from  this  failure  Is  that  the  LASA  m^'s  are  not  well  correlated 
with  spectral  shape  at  LASA.  We  did  not,  however,  design  our  experiment  to 
study  .these  problems  directly  and  do  not  yet  discount  completely  the  gains 
to  be  made  In  discrimination  by  correcting  m^  for  attenuation  differences, 
along  with  perhaps  regional  corrections.  These  corrections  will  require 
detailed  source-path  analysis  though. 

The  best  check  on  spectral  ratio  as  a discriminant  Involves  comparing 
explosions  with  earthquakes  In  close  enough  proximity  that  the  effect  due  to 
Q along  the  path  Is  essentially  unchanged.  We  have  examined  a small  body  of 
data  which  probably  fits  this  criterion,  namely  epicenters  In  the  southwest 
USSR,  and  found  that  the  separation  of  source  type  on  the  basis  of  the  LL 
spectral  ratio  Is  unsatisfactory.  This  result  Is  consonant  with  that  for 
earthquakes  and  explosions  In  the  NTS  area  and  further  degrades  the  standing 
of  P-wave  spectral  ratio  as  a reliable  discriminant  by  Itself.  However,  we 
have  seen  that  the  explosion  ratios  are  among  the  lowest  of  the  earthquake 
ratios;  and  on  this  basis  the  use  of  spectral  ratio  In  a multidiscriminant 
procedure  would  be  worthwhile  If  we  could  be  certain  that  the  dlscrlmlnantlon 
were  not  due  to  the  explosions  being  still  In  a source  region  of  higher  Q. 
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